Astronomy & Astrophysics manuscript no. aa'ms9922 


© ESO 2008 


September 11, 2008 





Photometry and polarimetry of the nucleus of comet 2P/Encke ^ 

bo"' 

O ' H. Boehnhardt\ G.P. Tozzi^, S. Bagnulo^, K. Muinonen"*, A. Nathues^ and L. Kolokolova^ 

o ■ 

CM . 

Oh 2 

D 3 
^ 4 



Max-Planck-Institute for Solar System Research, Max-Planck-Str. 2, D-37191 Katlenburg-Lindau, Germany (boehnhardt@mps.mpg.de) 
INAF Arcetri Observatory, Florence, Italy 

Armagh Observatory, College Hill, Armagh BT61 9DG, Northern Ireland, U.K. 
University Observatory, Helsinki, Finland 
University of Maryland, College Park, Maryland, USA 



' Received today / Accepted tomorrow 

^ : 

Cic ABSTRACT 



o 



Context. Short-period comet 2P/Encke is a favorable target for studies of light scattering by bare cometary nuclei. These studies enable 
assessment of the nucleus size, the material albedo, and the light-scattering properties of the cometary surface. 



Aims. We characterize the activity of the cometary coma and the physical properties of the nucleus of 2P/Encke such as its size, albedo, colors, 
, spectral slope, surface roughness, porosity, and single-particle properties. 

■ Methods. Broadband imaging photometry and broadband and narrowband linear polarimetry is measured for the nucleus of 2P/Encke over 
the phase-angle range 4-28 deg. An analysis of the point spread function of the comet reveals only weak coma activity, corresponding to 
a dust production of the order of 0.05 kg/s. The nucleus displays a color independent photometric phase function of almost linear slope (fi 
= 0.050±0.004 mag/deg). The absolute R filter magnitude at zero phase angle is 15.05±0.05, and corresponds to an equivalent radius for the 
nucleus of 2.43±0.06km (adopted albedo of 0.047). The nucleus color V-R is 0.47±0.07, suggesting a spectral slope S ' of 1 1 +8 %/lOOnm. The 
phase function of linear polarimetry in the V and R filter shows a widely color independent linear increase with phase angle (0. 12 ± 0.02 %/deg). 
We find discrepancies in the photometric and polarimetric parameters between 2P/Encke and other minor bodies in the solar system, which 
may indicate significant differences in the surface material properties and light-scattering behavior of the bodies. 

Results. The linear polarimetric phase function of 2P/Encke presented here is the first ever measured for a cometary nucleus, and its analysis 
encourages future studies of cometary nuclei in order to characterize the light-scattering behavior of comets on firm empirical grounds and 
provide suitable input to a comprehensive modeling of the light scattering by cometary surfaces. 
Conclusions. 
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1. Introduction 

After its discovery in 1786, the short-period comet 2P/Encke 
was observed during almost 60 perihelion passages and is one 
of the most significantly (more than 40 papers in refereed jour- 
nals) and longest (over 200 years) studied comets. With an or- 
bital revolution period of about 3.3 years, it is closer to the Sun 
than typical Jupiter family comets (JFC). Levison et al. (2006) 
suggested that 2P/Encke may be a rare example of a comet in 
an orbit detached from Jupiter's immediate control (like JFCs 
are) due to gravitational interaction with the terrestrial planets. 
They also argued that the nucleus may have been in a quies- 
cent state for a long part of the transition period into its current 
orbit. 

A series of scientific papers dealt with the nucleus proper- 
ties of 2P/Encke. Fernandez et al. (2000) measured an effective 



* Based on observations performed under program 078.C-0509 at 
the European Southern Observatory, Cerro Paranal, Chile 



nuclear radius of 2.4±0.3km, an axial ratio of at least 2.6, a ge- 
ometric albedo of 0.047 + 0.023, which is in the typical range 
for comets, and a rather steep photometric phase function with 
a linear slope of 0.06 mag/deg which was later confirmed by 
Jewitt (2002). A wider range of values for the effective radius 
(2.4-3.7 km) was determined from radar echoes of the comet 
(Harmon and Nolan 2005). Colors of V-R = 0.39+0.06 and B-V 
- 0.73+0.06 (Lowry and Weissman 2007) indicated only mod- 
erate spectral slope S' of the nuclear surface as found by Jewitt 
(2002) to be 8.9+1.6 %/lOOnm. The rotation period was a mat- 
ter of debate in various papers without a final conclusion (see 
Belton et al. 2005 and references therein). Attempts to deter- 
mine the orientation of the rotation axis were made by Sekanina 
(1991) and Festou and Barale (2000). Jewitt (2004) measured 
a slightly positive linear polarization of the nucleus at a single 
phase angle (polarization of 1.0-1.9 +0.5 % at 22°) in the vis- 
ible. The activity profile of the nucleus appears to be unusual 
with a coma being visible around aphelion (Fernandez et al. 
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2005), but no activity during its inbound arc as close as 1 .4 AU 
solar distance (Jewitt 2004). 

The likely absence of activity combined with the predicted 
nucleus brightness enable us to attempt successfully (i.e. with 
sufficient signal-to-noise and without significant coma con- 
tamination) polarization measurements of 2P/Encke. We re- 
port below on the results of the phase-angle-resolved linear 
polarimetry and quasi-simultaneous broadband photometry of 
this cometary nucleus in the visible wavelength range. The 
main aim of these measurements is to constrain the light- 
scattering properties of the nuclear surface, i.e. the single- 
scattering albedo and mean free path as a proxy for the typical 
grain size to be obtained through modeling. Furthermore, we 
could test the albedo-polarization relationship known from as- 
teroids (Zellner and Gradie 1976, Zellner et al. 1977, Lupishko 
and Mohamed 1996, Cellino et al. 1999), for a cometary nu- 
cleus. Our results are the first ever linear polarimetric observa- 
tions of a cometary nucleus over a wider phase-angle range. 

2. Observations and Data Reduction 

The observations of comet 2P/Encke were performed at Unit 
Telescope 2 of the Very Large Telescope observatory (VLT) 
at Cerro Paranal in Chile. The FORSl instrument was used 
in imaging and polarimetric modes with broadband Bessell V 
and R and two narrowband filters of central wavelengths and 
widths 485nm/37nm and 834nm/48nm, respectively; the filters 
are abbreviated below as "^^cont" for the blue and 'Vcont" for the 
red one. The narrowband filters were selected from the avail- 
able set of FORS filters such that their transmission ranges con- 
tained little to no contamination from emission bands of possi- 
ble coma gases, i.e. the narrowband filters cover mostly surface 
and dust reflected sunlight. The fecont and rcont filters were only 
used during the polarimetric measurements at the beginning of 
the program, since, due to the absence of significant coma con- 
tamination, we decided to maximize the signal-to-noise level 
by exclusive usage of the broadband V and R filters from 28 
Oct. 2006 onwards (see Table[T]l. 

2.1. Observations 

Service mode observations were the most appropriate choice 
given the specific requirements of this program: The short du- 
ration (1-1. 5h) of individual observing runs; a frequent, but ir- 
regular run schedule to cover the orbit of 2P/Encke at about 
regular phase-angle differences between 4 and 30 deg phase 
angle; and acceptable sky conditions for the observations. To 
the advantage of a dense phase-angle coverage, we relaxed the 
requirements for seeing (2") and sky transparency (thin cutus 
or better; in most cases the measurements were performed with 
clear sky) for our program. The comet was positioned at the 
center of the instrument field of view, a location that placed it 
automatically in the center of the central strip of the FORSl 
polarimetric mask used to separate the two beams of polarized 
light produced by the Wollaston prism in polarimetric mode. 
Exposure series of linear polarimetry of the comet were taken 



for each filter at A/2 retarder plate settings of position angles 
(with respect to North Celestial Meridian) 0°, 22.5°, 45°, 67.5°, 
90°, 112.5°, 135°, and 157.5°. The filter imaging and polari- 
metric exposure series were obtained subsequently within a 
few minutes and less than 1 hour, respectively, per observing 
run. Differential autoguiding on a nearby star at the speed of 
the comet in the sky was applied during the observations and 
the usual calibration exposures according to the FORS instru- 
ment calibration plan were taken (i.e. for imaging: photometric 
standard star field, bias and sky flatfield exposures; for the po- 
larimetry: bias and screen flatfield exposures plus a polarized 
and unpolarized standard star). The observing program was ex- 
ecuted in 9 runs during 1 Oct. and 14 Dec. 2006. Details of the 
observations are summarized in Table [1] 



2.2. Reduction of tlie Photometry Data 

All images were bias subtracted, and images used for photom- 
etry were then divided by a master flat field obtained from four 
sky flats taken at twilight. For the background subtraction, we 
subtracted in a first step a constant value that was estimated in 
regions far from the comet photometric center, where the con- 
tribution of a possible coma was negligible. A further correc- 
tion was evaluated and subtracted by measuring the EAf func- 
tion, defined by Tozzi et al. (2004), versus the projected nucle- 
ocentric distance, p. The SAf function is defined by the product 
of the geometric albedo multiplied with the total area covered 
by the solid component (usually the cometary dust; A'Hearn et 
al. 1995) in an annulus of radius p and unitary depth (Tozzi et 
al. 2004). For a normal comet, the EAf function is constant, and 
in the case of no coma it should be zero. Following a trial and 
error procedure, different background values were subtracted 
until the EAf function was measured to be constant versus p. 
As shown in Figure 2, EAf shows a peak at the cometary center, 
where the light-scattering contribution from the nucleus domi- 
nates, and decreases to values close to zero at larger projected 
distances from the nucleus (for a more detailed discussion see 
below). Finally, the images were calibrated both in magnitude 
and in Af. The value Af describes the product of the mean 
albedo A and the filling factor f of the dust grains in a measure- 
ment aperture of radius p. The product Afp is a measure of the 
dust production of the comet (for details on Afp see A'Hearn et 
al. 1995). Using standard star images we determed photometric 
zeropoints for the respective observing nights assuming atmo- 
spheric extinction coefficients of 0.1 14 and 0.065 mag/airmass 
and instrument color coefficients of 0.03 and 0.06 (plus solar 
colors for the comet to a first approximation) for V and R filters, 
respectively. The final two photometric parameters could not be 
determined from the available calibration images directly and 
were adopted after a careful analysis of the respective informa- 
tion on ESO's data quality information web page for the FORS 
instrument The cometary photometry was measured using a 
constant aperture diameter of 6". 



' see http://www.eso.org/instruments/l for technical details 



- see http://www.eso.org/observing/dfo/quaIityl 
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Table 1. Observing log and measurement results of the photometry and linear polarimetry of comet 2P/Encke. 



EPOCH 


PHASE 


r 


A 


Filter 


Ho 




Pu 


Pl 




Notes 


(UT) 




(All) 


(All) 




(mag) 


) 


('/() 


( '/< ) 


(') 




2006-10-01 01:40:02 


3.75 


2.75 


1.76 


V 




-0.90 + 0.51 


0.03 + 0.51 


0.90 


89+17 




2006-10-02 01:48:31 


4.01 


2.74 


1.75 


V 




-0.70 + 0.58 


-1.20 + 0.59 


1.39 


120 + 17 




2006-10-02 02:04:20 


4.02 


2.74 


1.75 






-0.78 + 0.42 


-0.42 + 0.43 


0.89 


104 + 19 




2006-10-02 02:20:12 


4.02 


2.74 


1.75 


R 




-0.68 + 0.32 


-0.07 + 0.32 


0.68 


93 + 15 




2006-10-02 02:40:03 


4.02 


2.74 


1.75 


^coat 




-1.45 + 0.40 


-1.24 + 0.40 


1.91 


110 + 9 




2006-10-02 02:56:20 


4.03 


2.74 


1.75 


R 


15.16 + 0.07 












2006-10-02 02:57:55 


4.03 


2.74 


1.75 


V 


15.66 ± 0.07 












2006-10-12 03:47:04 


8.09 


2.66 


1.71 


V 




-0.28 + 0.36 


-0.20 + 0.35 


0.35 


107 + 41 




2006-10-12 04:03:19 


8.10 


2.66 


1.71 


^cont 




-0.50 + 0.30 


0.13 + 0.45 


0.52 


83 + 29 




2006-10-12 04:19:47 


8.10 


2.66 


1.71 


R 




-1.12 + 0.32 


-0.52 + 0.33 


1.23 


102 + 10 




2006-10-12 04:40:12 


8.11 


2.66 


1.71 


^ cont 




-0.95 + 0.49 


0.02 + 0.98 


0.95 


89 + 30 




2006-10-12 04:56:35 


8.12 


2.66 


1.71 


R 


15.45 ± 0.07 












2006-10-12 04:58:06 


8.12 


2.66 


1.71 


V 


15.93 ± 0.09 












2006-10-12 05:02:34 


8.12 


2.66 


1.71 


R 


15.45 + 0.07 










[1] 


2006-10-12 05:04:11 


8.12 


2.66 


1.71 


V 


15.98 + 0.08 










[1] 


2006-10-20 00:21:07 


11.79 


2.60 


1.71 


R 


15.68 ± 0.07 












2006-10-20 00:22:28 


11.79 


2.60 


1.71 


V 


16.12 + 0.07 










[2] 


2006-10-20 00:33:30 


11.79 


2.60 


1.71 


V 




0.07 + 0.42 


-0.36 + 0.42 


0.37 


141 + 38 




2006-10-20 00:49:18 


11.80 


2.60 


1.71 


^cont 




-0.46 + 0.71 


0.43 + 0.66 


0.63 


68 + 44 




ZUUo-iU-zU Ul.Uj.io 


1 1 on 


Z.oU 


1. / 1 


K 




— U.Z4 + U.jZ 


U.Uj + U.ji 


U.Z4 


oj + 4J 




2006-10-20 01:25:05 


11.81 


2.60 


1.71 






-1.54 + 0.82 


-0.01 + 0.67 


1.54 


90+ 12 




2006-10-28 02:43:33 


15.46 


2.53 


1.75 


R 




-0.31 + 0.15 


0.08 + 0.15 


0.32 


83+17 




2006-10-28 03:18:40 


15.47 


2.53 


1.75 


V 




0.37 + 0.20 


-0.30 + 0.19 


0.48 


160 + 17 




ZUUo-lU-Zo Uj.j/.U/ 


1 j.4o 


Z.ji 


L./J 


K 


ij.oj + U.U/ 












2006-10-28 03:38:38 


15.48 


2.53 


1.75 


V 


16.17 + 0.07 












2006-11-04 00:41:39 


18.32 


2.47 


1.73 


R 




0.92 + 0.33 


-0.13+0.34 


0.93 


176 + 12 




2006-11-04 00:59:12 


18.32 


2.47 


1.73 


V 




0.89 + 0.43 


-0.30 + 0.43 


0.94 


171 + 17 




2006-11-04 01:28:12 


18.33 


2.47 


1.73 


R 


15.97 + 0.14 












2006-11-04 01:29:45 


18.33 


2.47 


1.73 


V 


16.40 + 0.22 












2006-11-07 00:39:12 


19.47 


2.44 


1.74 


R 




0.56 + 0.21 


-0.14 + 0.19 


0.58 


173 + 12 




2006-11-07 01:14:18 


19.47 


2.44 


1.74 


V 




1.13 + 0.25 


-0.12 + 0.25 


1.13 


177 + 7 




2006-11-07 01:34:45 


19.47 


2.44 


1.74 


R 


15.97 + 0.10 












2006-11-07 01:36:19 


19.47 


2.44 


1.74 


V 


16.58 + 0.10 












2006-11-10 00:56:40 


20.54 


2.42 


1.75 


R 




1.14 + 0.22 


-0.02 + 0.22 


1.14 


179 + 6 




2006-11-10 01:32:37 


20.55 


2.42 


1.75 


V 




0.55 + 0.25 


0.32 + 0.23 


0.64 


15+15 




2006-11-10 01:53:11 


20.55 


2.42 


1.75 


R 


15.98 + 0.08 












2006-12-14 01:08:21 


27.97 


2.09 


1.94 


R 




1.90 + 0.40 


-0.37 + 0.41 


1.93 


175 + 7 


[3] 


2006-12-14 01:15:21 


27.97 


2.09 


1.94 


V 




1.20 + 0.96 


-1.51 + 1.12 


1.93 


154 + 22 


[3] 



Explanations: The table lists the observing epoch (midpoint of exposure series), the Sun (r) and Earth (A) distances, the phase angle (0), the 
filter used, the absolute filter brightness m( 1 , 1 ,</>) for r = A = 1 AU and phase angle (p, the Stokes parameters Pq and Pu , the fraction of linear 
polarization /"l and the angle ^ of maximum polarization. Comments are indexed and explained below: [1] Photometric exposures in V and R 
were repeated because of guide probe vignetting. [2] No zeropoints available in V on 2006-10-20. We used those obtained on 2006-10-19. [3] 
Target very close to a bright star, accurate polarimetry and photometry impossible. 



2.3. Reduction of the Linear Polarimetry Data 



All polarimetric images were bias subtracted and flatfield cor- 
rected. For the background subtraction, we followed the same 
strategy used for the imaging data of the comet, by treating the 
ordinary and extra-ordinary beams separately. A further refine- 



ment in the background estimate was obtained by measuring 
(for each position angle of the retarder waveplate) the ratio 



(1) 



where f° is the flux in the ordinary beam and the flux in 
the extraordinary beam, in annuli of increasing radius. We then 
readjusted the value used for background subtraction until the 
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measured h values became independent of the projected nucle- 
ocentric distance p. In the R band, typically, the background 
values obtained after the first step (i.e. estimated from regions 
far from the cometary nucleus) were approximately 800 e ; the 
first correction, based on the analysis of the EAf function, was 
approximately 10 - 20 e ; the final correction derived by in- 
sisting that the ratio h was constant with p, was approximately 
1 - 2e-. 

Linear polarization was then calculated from 
^' = -2r(a,)cos(4a,), (2) 

1=0 

and 

2 " 

m' = -2r(a,)sin(4a,), (3) 

1=0 

where q' and u' are the reduced Stokes parameters defined ac- 
cording to Shurcliff (1962) and measured using as the refer- 
ence direction the North Celestial Meridian, and a,- = 22.5° X ; 
are the position angles of the retarder waveplate. From q' and 
u', we finally calculated the linear polarization Pq = Q/I 
and Pu = U/I, by adopting as the reference direction the 
perpendicular to the great circle passing through the comet 
and the Sun at the observing epoch, as explained in Landi 
Degl'Innocenti et al. (2007). Pq and Pu allow one to deter- 
mine the total fraction of linear polarization f l - -^jP^g + 
and the angle of maximum polarization f . Pq represents the 
flux perpendicular to the plane Sun-Comet-Earth (the scatter- 
ing plane) minus the flux parallel to that plane, divided by the 
sum of the two fluxes, and ^ is the angle between the perpen- 
dicular to the scattering plane and the direction of maximum 
polarization. As for the photometry, a fixed aperture of 15" di- 
ameter was applied for the polarimetric measurements of the 
comet. In the presence of coma, we have implicitly assumed 
that the polarization of the coma does not depend on p (see 
below for a discussion of the influence of a faint coma on the 
polarization data of the comet). 

3. Results 

Search for coma: No coma is detectable by direct visual in- 
spection of the comet images. To confirm the presence of a faint 
coma possibly undetectable by visual inspection and measure 
its contribution to the object brightness, two numerical analysis 
methods were applied. The first and classical method (see for 
instance Boehnhardt et al. 2002) compared the azimuthally av- 
eraged profiles of the comet with the instrumental point spread 
function (PSF) measured from background stars in the same 
exposure. Since the comet had a proper motion and the tele- 
scope was tracking the comet, the stars appeared as short trails. 
Hence, the stellar PSF can only be measured in one direction 
i.e. perpendicular to the cometary motion and averaged along 
the trail direction. Figure[T]shows that the stellar and cometary 
PSF are almost identical, in particular in the PSF wings and 
close to the background level at which the coma light is ex- 
pected to provide the most significant contribution. 
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Fig. 1. Azimuthally averaged PSF profile of the comet (dashed 
line) compared with that of a background star in the same ex- 
posure. The image was taken on 28 Oct. 2006 through the V 
filter. 

To increase the sensitivity to coma detection, all polarimet- 
ric images recorded in a single filter during a single observ- 
ing night were centered on the comet position and coadded, 
reaching a total integration time of approximately 1000s (in- 
stead of the short 60s exposures used for normal photomet- 
ric observations). Due to seeing variations during exposure se- 
ries, direct comparison of the PSFs of comet and background 
stars was complicated and usually impossible. For these com- 
posite images a second analysis method was therefore applied: 
Measurement of the EAf as function of p (see above). In the 
case of the absence of coma, this function should approach zero 
at radial distance p of about 3 FWHM from the central bright- 
ness peak, with FWHM being the value of the full width at half 
maximum of the PSF of the coadded comet image. In the case 
of the coadded polarimetric images of 2P/Encke, all EAf pro- 
files show small, but significant non-zero fluxes at larger dis- 
tances from the center of the coadded comet image. We recall 
that the background level has been subtracted using EAf as de- 
scribed above and that any residual background level would re- 
sult in a (likely linear) non-constant trend of SAf versus p. The 
fact that for the 2P/Encke polarimetric images SAf is constant 
with p means that the intensity profile of the comet declines 
with 1/p, a phenomenon that is typical of an expanding dust 
coma around the nucleus. We therefore conclude that a weak 
coma was present around the comet during the entire period of 
our observations, i.e. from early Oct. to mid Dec. 2006 when 
the comet was approaching the Sun from 2.75 to 2. 1 AU. Given 
the overall SAf profile, it is unlikely that the central brightness 
peak of the comet images is due to an unresolved dense dust 
coma; instead, we consider the photometric flux of the peak 
produced by surface-reflected sunlight. 

Flux calibration of the polarimetric images was obtained 
by comparing the inner part of the SAf profile of the coadded 
polarimetric image with that of a calibrated normal filter ex- 
posure of the comet taken in the same filter during the same 
night. Figure |2] shows an example of the EAf profile of a coad- 
ded polarimetric image series of comet 2P/Encke. Within the 
error margins, the ZAf profiles obtained for the various observ- 
ing nights show no trend, neither with heliocentric distance 
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nor with phase angle. The average level of the weak coma 
flux corresponds to Afp values (for the definition of Afp see 
above) of 0.65+0.35 cm in V and 0.49+0. 19 cm in . They are 
equivalent to a dust production rate Qdust of the order of about 
0.05 kg/s, assuming a simple empirical relationship between 
Afp and Qdusi used by Kidger (2004). The measured Afp val- 
ues equal no more than 1 percent of the total cometary signal 
and its effect on the photometric and polarimetric data analysis 
of 2P/Encke described below, can be neglected. We also note 
that during the reduction of the polarization images a weak ex- 
tended coma was subtracted as part of the general background 
signal measured beyond 3 FWHM of the PSF of the comet. 
Only second order contributions to the Q and U measurements 
of the comet signal may remain from polarized light of the dust 
coma around the nucleus. Since estimated to be below 1 per- 
cent of the measured Stokes parameters, we ignored these in the 
overall error analysis of the polarimetric results of 2P/Encke. 
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Fig. 2. SAf profile as function of the nucleocentric distance 
p, obtained from polarimetric observations in V filter on 28 
Oct. 2008. The total exposure time is 960s. For radial distances 
greater than 2000 km, ZAf is small but different from the back- 
ground level of zero flux. Hence a very faint coma is present 
around the cometary nucleus. 

Photometric phase function of the nucleus: Figure [3] 
shows the photometric phase function of the nucleus of comet 
2P/Encke for broadband V and R filters. The values plotted 
therein, the so-called m(l,l,(^) magnitudes (see Table [T]i, are 
derived from the observed filter magnitudes by removing the 
dependencies on the Sun and Earth distances for the respec- 
tive observing epochs. The figure shows a steep almost linear 
brightness increase for the nucleus with decreasing phase an- 
gle. It also suggests that there is - at least in the R filter data 
- a slight non-linear increase in brightness toward zero phase, 
starting all the way from = 20°. Although this deviation from 
linear slope is within the expected amplitude of brightness vari- 
ations due to the rotation of the non-spherical nucleus, its sys- 
tematics suggest that the opposition effect may play a role in the 
nucleus phase function of 2P/Encke. Ignoring the possible mi- 
nor non-linearity in the phase functions of the comet, we found 
almost identical slopes fi when fitting a linear phase curve, i.e. 
for the V filter data jS = 0.051+0.004 mag/°and for the R fil- 
ter ones p - 0.049+0.004 mag/°. Our p is slightly smaUer, 



although consistent within the eiTor margins, than the linear 
slope parameter {fi - 0.06 mag/°) determined by Fernandez 
et al. (2000) from a compilation of published nucleus magni- 
tudes of the comet measured over a wider phase-angle range 
and that found by Jewitt (2004) for 2P/Encke (fi = 0.060+0.005 
mag/°) derived from a much sparser dataset. We conclude that 
the phase function of 2P/Encke follows a color-neutral almost 
linear brightening law with /3 = 0.050+0.004 mag/°, which is 
slightly steeper than the canonical value (J3 = 0.04mag/°; Lamy 
et al. 2004), although within the range found for other cometary 
nuclei (Jewitt 2004, Lamy et al. 2004). 
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Fig. 3. The brightness phase function of the nucleus of comet 
2P/Encke in broadband V and R filters. The plot shows the mea- 
sured m( 1,1,0) brightness, in addition to the linear fit functions 
wersus phase angle. Filled squares = R filter, open triangles - 
V filter. 

Nucleus size and colors: Extrapolation of the linear phase 
function of 2P/Encke to zero phase angle provides the abso- 
lute magnitude m(l,l,0) of the nucleus, i.e. 15.50+0.06 and 
15.03+0.05 mag for V and R filter, respectively. With the ge- 
ometric albedo of 0.047 (Lamy et al. 2004; Fernandez et al. 
2000), we obtain an equivalent radius of 2. 43+0. 06km in R, in 
good agreement with the radius estimate by Fernandez et al. 
(2000; 2.4+0.3km), Kelley et al. (2006), and Campins et al. 
(1988) and slightly smaller than radar echo results by Harmon 
and Nolan (2005; 2.42 - 3.72km). We emphasize that our re- 
sults may be affected by insufficient sampling of the rotation 
light curve and may deviate from earlier findings because of 
different viewing aspects of the nucleus along the orbit and with 
time (see Belton et al. 2005). The V-R color of the nucleus is 
0.47+0.07 mag, which coiTesponds to an intrinsic color (i.e. 
corrected for solar V-R color) of 0. 1 1 +0.07 mag or a spectral 
slope S' of 1 1±8 %/lOOnm. Our mean V-R color and the spec- 
tral gradient S' are slightly higher than that measured by Jewitt 
(2004) and Lowry & Weissman (2007), although still within 
the error margins. 

Polarimetric phase function of the nucleus: Stokes Pq 
and Pu of the Hnear polarimetry measured for 2P/Encke are 
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compiled in Table [T] The Pq and P^ values listed refer to 
the scattering plane Sun-comet-observer. The polarization 
with respect to the scattering plane (listed in Table [TJ is plot- 
ted versus phase angle in Figure |4] for the two broadband fil- 
ters used (V and R). In both broadband filters the polarization 
shows the same, linear increase with increasing phase angle (p 
(from about -1% at 0=4°to about +1.1% at 0=28°). The min- 
imum polarization values for V and R correspond to the up- 
per limits only, since our polarization phase functions do not 
show the expected turn-over towards zero polarization for small 
phase angles (the minimum polarization is rather the lower end 
of the linear phase function slope and the turn-over occurs at 
(f> < 4°only). Zero polarization is passed at about 13°phase an- 
gle (inversion angle). The slope of the polarization increase is 
determined by linear regression to be 0.123 + 0.013 %/°in V 
and 0.129+0.021 %/°mR,i.e. within the errors the phase-angle 
gradients are identical and do not change with filter. 

From the values listed in Table[Tl it is obvious that for small 
Pq, i.e. close to the overall measurement uncertainty. Stokes 
Pu has comparable amplitudes to Pq which produces a con- 
siderable scatter and large uncertainties in the values for the 
position angle ^ of maximum linear polarization. The position 
angle ( is well determined for higher polarization values and 
agrees well with the position angle of the light scattering plane. 
However, the close-to-zero level of Stokes Pu for all measure- 
ments provides evidence for the correctness of the data reduc- 
tion, and is in agreement with theoretical expectations for a 
light-scattering-plane oriented polarization (see Muinonen et 
al. 2002). For the narrowband measurements using filters fecont 
and Tcont, Stokes Pq is systematically higher for the red con- 
tinuum filter compared with the blue, although all results are 
affected by relatively large errors. Within the margins of the 
errors the measured Stokes parameters agree with the Stokes 
parameters measured for the V and R filters, respectively. 



3.0 - 
2.5 

2.0 - 




-2.0 -\ ^ > > > > 1 

5 10 15 20 25 30 
Phase Angle phi (deg) 

Fig. 4. Linear polarization P of comet 2P/Encke versus phase 
angle. Filled triangles/squares = V/R filter measurements, 
long/short dash lines = Unear fits to V/R filter measurements. 



To our knowledge, 2P/Encke is the only comet that has lin- 
ear polarization at the nucleus measured. Jewitt (2004) pub- 
lished a set of polarization values (1.5, 1.0, 1.9% ±0.5 %), ob- 
tained on 22 Oct. 2003 at phase angle 22.1 "through cometary 
dust continuum filters (central wavelengths of 445.3, 525.9, 
713.3 nm). His results are consistent with our own. 

Test of the asteroid polarization-albedo relationship: 
We also applied the polarization-albedo relationship for as- 
teroids (Zellner and Gradie 1976, Zellner et al. 1977) to the 
2P/Encke polarimetry data to validate this empirical method 
of albedo estimation for the nucleus of the comet. Using the 
phase function parameters of Lupishko and Mohamed (1996) 
we derive a geometric albedo of 0.145 (range 0.13-0.16), when 
using the average slope of the polarization phase function (see 
above), and of 0.08 (range 0.06-0.16) when using the values 
of minimum polarization P,„,„ - which by itself is an upper 
limit only. Unfortunately, in our measurements of 2P/Encke 
P„„„ is not well determined and the upper limit used instead 
shows relatively large errors (see Table [TJ. In conclusion, the 
empirical polarization-albedo relationships for asteroids sug- 
gests a significantly higher geometric albedo for the nucleus 
of 2P/Encke compared to that obtained from the more realistic 
approach of visible and thermal IR measurements of the nu- 
cleus (0.047; Lamy et al. 2004; Fernandez et al. 2000). This 
may indicate that either the empirical polarization-albedo re- 
lationship of asteroids does not apply to cometary nuclei, or at 
least not for 2P/Encke, or, if applicable, that it requires different 
fitting parameters than determined from asteroids. It is indeed 
been proposed that the empirircal polarization-albedo rule is 
not applicable to extremely dark objects such as cometary nu- 
clei (Dollfus & Zellner 1979). The negative outcome of this 
test may imply that the surface constitution of 2P/Encke differs 
from those of main-belt and near-Earth asteroids for which the 
empirical relationship is calibrated by measurements. 

4. Discussion and Conclusions 

We have presented photometric and linear polarimetric phase 
functions of 2P/Encke measured during the 2006 approach of 
the comet to the Sun. The weak coma found about the comet 
does not affect significantly the photometric and polarimetric 
light-scattering phase functions measured, i.e. the phase func- 
tions reflect the light-scattering behavior of the nucleus sur- 
face. Both phase functions show an almost linear behavior with 
phase angle <p (over the measured range from 4 to 28 deg). 
Trends in the phase function with wavelength are not obvious 
for the nucleus photometry, but may exist for the linear po- 
larization. A small opposition surge is tentatively suggested in 
the photometric data, but needs verification by further obser- 
vations at smaller phase angles. Similarly, the value of mini- 
mum polarization P„„„ in the polarization phase curve should 
be reassessed by new observations covering the range to zero 
phase angle that is of immediate interest for modeling the op- 
position effects in the light-scattering by the cometary nucleus. 
Our radius estimation and spectral slope S' agree with earlier 
measurements, in particularly well with those obtained from IR 
observations. The low activity level of the nucleus at solar dis- 
tances of 2.75 to 2.1AU suggests that the cometary nucleus is 
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Table 2. Comparison of photometric and polarimetric properties of comet Encke and other small bodies. 



Object 


Albedo 


Phase coefficient. 


Spectral 






Pol. slope 


Pol.Inv. 


Pol.spect. 




Pv 


/? for 5 < < 30° 


Slope 


(%) 


(°) 


(%/°) 


angle 0„ 


gradient 






(mag/°) 


(%/lOOnm) 








n 


at < 30" 


2P/Encke 


0.047 


0.048 


11 


> 1.4 


<4 


0.123 


13 


negative 


















decreases with (p 


Comet nuclei 


0.02-0.06 


0.043-0.059 


0-41 


7 


7 


7 


7 


7 


Comet dust 


0.05 


0.018-0.035 


-1-12 


1.5 


9.6 


0.31 


21.8 


negative 


















decreases with (f) 


2060/Chiron 


0.07 


0.045 


1 


1.5 


1.5-2 


Out of range 


9(?) 


negative ? 


C asteroids 


0.02-0.09 


0.042 


-3-20 


1.73 


9.4 


0.28 


20.5 


positive 


















increases with (p 


S asteroids 


0.085-0.33 


0.029 


-15-88 


0.77 


6.3 


0.09 


20.1 


negative 


















increases with (p 


E asteroids 


0.34-0.6 


0.019 


rather flat spectra 


0.3 


5.8 


0.04 


17.8 




M asteroids 


0.074-0.25 


0.033 


6-15 


0.96 


7.7 


0.09 


23.5 


negative 


















increases with if> 


F asteroids 


0.024-0.086 


0.044 


-24 (Interamnia) 


1-1.4 


7.6-8 


0.35-0.4 


14-17 


slightly negative or 








bluish spectra 










neutral 



Explanations: Column "Spectral Slope" lists the slope of the visible spectrum of the objects as determined from measured filter colors 
corrected for the color of the Sun (see Boehnhardt et al. 2001); "l/'mml" gives the value of minimum polarization, "^m,„" the phase angle 
of \Pmw\, "Pol. slope" the gradient of the polarimetric phase function and "Pol.Inv." the polarization inversion angle, i.e. when the sign of 
polarization Pt changes from negative to positive. Column "Pol.spect. gradient at (p < 30°" describes the overall spectral shape of the object 
spectra for phase angles below 30 °. The references for the table values are: albedo for comet nuclei - Lamy et al 2004; filter colors of objects 
- Snodgrass et al. 2008; phase coefficients from space missions - Soderblom et al. 2002, Li et al. 2006 phase coefficients from ground-based 
observations - Jewitt and Sheppard 2004, Jewitt et al. 2003; albedo for 2060/Chiron - Groussin et al. 2004; V - R color of 2060/Chiron - Davies 
et al. 1998; the rest of 2060/Chiron's data - Bagnulo et al. 2006; asteroid albedos - Tedesco et al. 1989; asteroid phase coefficients - Bowell 
et al. 1979; asteroid filter colors - Hansen 1976, Chapman & Gaffey 1979; polarimetric asteroid data - Lupishko 2006; spectral gradients of 
polarization for asteroids - Belskaya et al. 2008; comet dust data - Kolokolova et al. 2004, Kiselev & Rosenbush 2004, Meech & Jewitt 1987. 
Question marks indicate that no comparison data are available. 



covered by a crust that prevents major outgassing. This appears 
to be true during this part of the orbit, even though 2P/Encke 
has entered the distance range where water sublimation usually 
generates significant cometary activity. 

Table |2]compares light-scattering properties of the nucleus 
of comet 2P/Encke with those of other solar-system bodies, 
that may have similar surfaces and were observed in similar 
conditions (phase angle, wavelength). The majority of the data 
are for V filter - if not indicated otherwise in the table - at the 
phase angles 5-30°. The analysis of Table|2]indicates that comet 
2P/Encke displays photometric properties typical of cometary 
nuclei. As for other cometary nuclei, the slope of the phase 
function differs from that measured for cometary dust, which 
might be due to the diff'erence between light scattering on in- 
dividual dust particles and their compact layers at the nuclear 
surfaces. 

Among asteroids, the light scattering behavior of C-type 
objects resembles most closely that of comet 2P/Encke, at 
least its photometric properties. There are major differences 
between the polarization phase curve (Figure |4]i of the nu- 
cleus of 2P/Encke and the polarization phase curves of aster- 
oids of diff'erent taxonomy. Primitive asteroids of type C, P, 
G, and B and evolved asteroids of type V, S, M, E, and A 
have, in general, larger inversion angles than 2P/Encke (cf. 



Muinonen et al, 2002; Fornasier et al., 2006), while only primi- 
tive F-type asteroids exhibit similar inversion angles (Belskaya 
et al., 2005). However, other polarimetric characteristics of F- 
type asteroids differ significantly: The angle of the polariza- 
tion minimum is at least twice as large, and the slope of polar- 
ization at the inversion angle is almost three times larger than 
that of comet 2P/Encke. Moreover, F-type asteroids do not ex- 
hibit red colors, but display neutral or even bluish surface col- 
ors. Unfortunately, we cannot compare the properties of comet 
2P/Encke with icy bodies such as Transneptunian objects or the 
small satellites of Saturn, Uranus, and Neptune since these ob- 
jects can be observed only at small phase angles. The Centaur 
2060/Chiron may be the object most appropriate for the com- 
parison. Compared to comet 2P/Encke, it has however nei- 
ther similar photometric (too high albedo, too red color) nor 
polarimetric data (too small inversion angle). The polarime- 
try therefore reveals unique properties of the nucleus of comet 
2P/Encke that may be typical for cometary nuclei but differ sig- 
nificantly from the properties of other solar-system bodies. This 
uniqueness does not manifest itself photometrically. It may in- 
dicate a specific composition or structure of the surface layer 
of cometary nuclei. Whether the narrow negative polarization 
branch can be explained qualitatively by the single-scattering 
interference mechanism (Muinonen et al. 2007) remains to be 
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answered by a future study. Based on that mechanism, nega- 
tive polarization branches of icy objects which corresponds to 
smaller real parts of refractive indices, can be expected to be 
narrower than those of siUcate-rich stony objects, which repre- 
sents larger real parts of refractive indices. 

The polarization minimum of 2P/Encke, which is not well 
constrained by the present data set, is less than or equal to - 
1 % and, as mentioned above, it is possibly within the range of 
data for F-type asteroids that have polarization minima between 
about -1.0 and -1.5% (Belskaya et al. 2005). However, we reit- 
erate that the polarization phase curve of 2P/Encke is the only 
data set for comets available to ourselves. The number of F- 
type asteroids with measured polarization phase curves is also 
small such that, in terms of polarization properties, it is difficult 
to decide whether comets and F-type asteroids are similar. 

Our 2P/Encke results have shown that the empirical albedo- 
polarization relationships of asteroids cannot easily be applied 
to cometary nuclei, possibly because the surface constitution 
and light-scattering properties are different. Whether such a re- 
lationship can be established for cometary nuclei remains un- 
decided until more objects are studied and a self-consistent 
modeling of photometric, spectroscopic, and polarimetric light- 
scattering properties for comets is available. 
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